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I .  INTRODUCTION 


This  report  addresses  three  main  areas  of  concern  to  STORM  model 
users.  These  are: 

1.  Collection  of  site-specific  data  for  calibration 

2.  Model  calibration 

3.  Model  application 

The  recommendations  presented  here  are  based  on  the  writer's  experience 
In  applications  on  several  studies  conducted  at  the  Hydrologic  Engineering 
Center  (HEC).  The  studies  were  sponsored  by  several  Corps  of  Engineers 
district  and  Division  offices. 


II.  SUMMARY 

'This  report  provides  specific  Information  on  calibration  and  appli¬ 
cation  of  the  Storage,  Treatment,  Overflow  Runoff  Model  (STORM).  The  STORM 
model  is  Intended  for  use  In  simulation  of  the  quantity  and  quality  of 
storm  water  runoff.  In  particular,  the  report  discusses  procedures  for 
collection  of  rainfall,  runoff  quantity  and  quality  data.  Procedures  are 
recommended  for  management  of  the  collected  data.  Recommendations  are 
provided  for  use  of  the  site-specific  data  In  calibration  of  the  model. 

The  calibrated  model  can  then  be  used  for  two  Important  planning  components 
of  a  storm  water  study.  These  are: 

1.  Prediction  of  wet-weather  pollutographs  (mass  loading  curves) 
for  use  In  a  receiving  water  assessment  model.  These  polluto- 
graphs  can  Include  both  surface  runoff  and  dry  weather  flow  In 
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combined  systems.  Since  the  computations  are  based  on  land  use, 
the  Impact  of  land  use  change  can  be  evaluated. 

2.  Preliminary  sizing  of  storage  and  treatment  facilities  to 
satisfy  desired  criteria  for  control  of  storm  water  runoff. 

The  model  will  analyze  a  matrix  of  combinations  of  storage  and 
treatment  rates.  Results  Include  frequency  Information  on 
quantity  and  quality  of  washoff  of  pollutants  and  soil  erosion, 
as  well  as  frequency  Information  on  the  quantity  and  quality  of 
storage  overflows. 


III.  BACKGROUND  OF  THE  STORM  MODEL 


The  quantity  portion  of  STORM  was  developed  for  the  City  of  San 
Francisco  by  Water  Resources  Engineers,  Inc.,  (WRE)  of  Walnut  Creek, 

California.  The  water  quality  computations  were  added  In  1973  by  WRE 
while  under  contract  with  the  HEC.  Since  then,  the  HEC  has  added  other 
capabilities  including  snowmelt  and  land  surface  erosion  computations  and 
prespecified  hydrographs.  Resource  Analysis,  Inc.,  of  Cambridge,  Massachusetts, 
added  the  capability  to  simulate  the  quantity  and  quality  of  dry  weather  flow. 
The  HEC  added  the  U.S.  Soil  Conservation  Service  (SCS)  Curve  Number  Technique 
for  runoff  computation  and  the  SCS  unit  hydrograph  technique.  A  future 
version  will  Include  channel  routing  and  combining,  and  a  "planning  level" 
stream  water  quality  module. 
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IV.  STORM  CONCEPTS 


r 


1.  General .  STORM  Is  a  continuous  simulation  model  that  can  be  used 
for  prediction  of  the  quantity  and  quality  of  storm  Mater  and  dry  weather 
flow  (domestic,  commercial,  industrial  and  pipe  infiltration).  The  model 
provides  analyses  that  can  be  used  to  satisfy  two  primary  (but  related) 
study  objectives.  These  are  (1)  prediction  of  wet  weather  pollutographs 
that  can  be  used  In  a  receiving  water  assessment  model  and  (2)  provide 
statistical  Information  to  aid  In  the  selection  of  storage  capacities  and 
treatment  rates  required  to  achieve  desired  control  of  storm  water  runoff. 

Wet  weather  pollutographs  can  be  predicted  for  Individual  historical 
events.  The  pollutographs  consist  of  hourly  rainfall,  runoff,  masses  of 
pollutants,  and  pollutant  concentrations.  This  Information  can  be  used 
directly  by  a  receiving  water  assessment  model. 

Another  portion  of  the  output  from  STORM  provides  statistical 
Information  based  on  an  analysis  of  the  hourly  precipitation  record. 

Statistics,  such  as  average  annual  runoff,  average  annual  washoff  of  each 
pollutant,  average  annual  overflow  from  storage,  and  average  annual  pollutant 
overflow  from  storage,  are  provided.  This  Information  can  be  used  to  aid  the 
selection  of  storage  capacities  and  treatment  rates  required  to  achieve 
desired  control  of  storm  water  runoff. 

Runoff  quantity  can  be  computed  by  one  of  three  methods,  the 
coefficient  method,  the  SCS  curve  number  technique  or  a  combination  of  the 
two.  In  the  coefficient  method,  average  annual  runoff  coefficients  for  the 
pervious  and  impervious  areas  of  the  watershed  are  specified,  and  subsequently. 
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weighted  according  to  the  total  fraction  of  impervious  area  in  the  watershed 
so  as  to  obtain  a  single  composite  runoff  coefficient.  This  coefficient  Is 
then  used  for  each  rainfall  event  In  the  precipitation  record  to  calculate 
runoff  excess  above  depression  storage  regardless  of  rainfall  or  soil 
characteristics.  This  method  may  not  produce  accurate  or  properly  shaped 
hydrographs  for  individual  rain  events,  but,  when  calibrated,  may  produce 
sufficiently  accurate  volumes  of  runoff. 

The  SCS  Curve  Number  Technique  uses  a  simple  curve  to  relate 
accumulated  runoff  to  accumulated  rainfall.  The  curve  number  Is  related 
to  the  soil  type  and  antecedent  moisture  conditions.  The  procedure  Includes 
use  of  an  Initial  abstraction  (depression  storage)  variable  wMch  must  be 
exceeded  before  any  runoff  can  occur  for  a  given  storm.  Thereafter,  the 
program  operates  on  the  curve  for  determination  of  runoff.  The  curve 
approaches  a  45  degree  slope,  i.e.,  near  zero  Incremental  Infiltration 
would  occur  at  the  end  of  a  very  large  storm.  Since  STORM  Is  a  continuous 
simulation  model,  only  initial  curve  numbers  (for  each  land  use)  are 
required.  The  model  computes  the  soil  moisture  storage  capacity  at  the 
beginning  of  each  storm  In  the  record  based  on  recovery  of  soil  moisture 
capacity,  initial  abstraction  and  percolation  during  dry  periods.  The 
curve  number  is  expressed  in  terms  of  Inches  of  soil  moisture  storage  for 
input  to  the  model  and  for  computations. 

The  combination  method  uses  the  SCS  method  an  pervious  areas  and 
the  coefficient  method  on  Impervious  areas  of  the  watershed. 

Storm  water  quality  Is  computed  by  using  an  exponential  washoff 
equation.  The  equation  relates  the  mass  of  pollutants  washed  off  during 


each  hour  to  the  current  mass  of  pollutants  on  the  watershed,  the  runoff 
rate,  and  an  exponent  governing  the  rate  of  pollutant  washoff. 

There  are  two  methods  by  which  the  pollutant  accumulation  may  be 
simulated.  In  the  dust  and  dirt  method,  the  mass  of  each  pollutant  on  the 
watershed  Is  computed  as  a  fraction  of  the  net  accumulated  dust  and  dirt  on 
each  land  use  at  the  beginning  of  each  period  of  rainfall.  The  dust  and 
dirt  accumulates  at  a  linear  rate  for  each  land  use  during  dry  hours. 
Calibration  of  this  method  Involves  adjustment  of  the  dust  and  dirt 
accumulation  rates,  the  pollutant  fractions  of  the  dust  and  dirt,  and  the 
washoff  coefficient  so  that  the  predicted  pollutant  concentrations  most 
nearly  match  those  from  measured  data.  In  the  dally  pollutant  buildup 
method,  the  mass  of  each  pollutant  Is  computed  by  a  pollutant  accumulation 
rate  for  each  constituent  In  terms  of  pounds/acre/day  for  each  land  use. 

The  pollutants  build  up  linearly  during  dry  hours.  Calibration  of  this 
method  Involves  adjustment  of  the  dally  pollutant  accumulation  rates  and 
the  washoff  exponent. 

Tne  overall  model  operation  Involves  the  interaction  of  8  main 
processes.  These  are:  precipitation,  runoff,  pollutant  accumulation, 
pollutant  washoff,  dry  weather  flow,  storage,  treatment,  and  overflow. 

The  program  computes  runoff  from  rainfall  or  rainfall  plus  snowmelt  and 
the  associated  pollutant  washoff.  Runoff  In  excess  of  the  specified 
treatment  rate  is  diverted  Into  storage  for  subsequent  treatment.  Runoff 
in  excess  of  both  the  treatment  rate  and  storage  capacity  Is  considered 
overflow  and  is  diverted  directly  into  the  receiving  waters.  The  program 
provlues  considerable  information  on  the  frequency,  quantity  and  quality 
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of  overflows.  Figure  1  shows  a  schematic  representation  of  the  major 
processes  modeled  by  STORM. 

The  reader  Is  referred  to  the  STORM  user's  manual  (Reference  2)  for 
details  of  the  program  operation,  Input  data  description,  and  samples  of 
input  and  output. 

2.  Limitations.  STORM  has  several  limitations  associated  with  the 
manner  in  which  calculations  are  made  and  the  nature  of  certain  Input 
variables.  The  model  uses  lumped  parameter  hydrologic  techniques  for 
computing  precipitation  excess  and  defining  subbasin  runoff.  Empirical 
techniques  are  used  for  computing  pollutant  washoff  and  land  surface  erosion 
and  normally  require  calibration  using  site-specific  data  for  the  study 
watershed. 

Another  limitation  is  the  one  hour  computation  interval.  The  main 
reasons  that  STORM  uses  a  one  hour  computation  Interval  are  that  continuous 
hourly  precipitation  records  are  usually  available  on  magnetic  tape  for 
most  locations  and  the  one  hour  interval  Is  probably  the  lower  limit  in 
terms  of  manageability  of  data  for  a  continuous  analysis  and  computer 
processing  time.  The  one  hour  Interval  affects  the  lower  limit  of  size  of 
subbasins  that  may  be  modeled.  One  should  not  expect  STORM  to  produce 
properly  shaped  hydrographs  for  subbasins  having  times  of  concentrations 
less  than  one  hour.  In  most  planning  studies,  however,  the  one  hour 
interval  should  not  represent  a  limitation  since  the  subbasin  sizes  are 
normally  large  enough  to  have  times  of  concentrations  greater  than  one  hour. 
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■Igurel.  MAJOR  PROCESSES  MODELLED  BY  STORM 


V.  INPUT  DATA  REQUIREMENTS 


1.  Summary.  The  following  Is  a  summary  of  Input  data  required  by 
the  major  program  options.  Chapter  VI  presents  detailed  Information  on 
collection  and  management  of  the  required  data  for  calibration. 

Data  Source 

1 .  Hydrological  Data 

Hourly  precipitation  U.S.  National  Weather  Service 

Evaporation  U.S.  National  Weather  Service 


a.  Coefficient  Method 

Depression  Storage  Hydrologists 

Runoff  coefficient  for 

pervious  areas  _  Analysis  of  historical  rainfall 

runoff  records 

Runoff  coefficient  for 
impervious  areas  _ 


b.  SCS  Method 


Maximum  soil  moisture 
retention  capacity 

Maximum  Initial  abstraction 
capacl ty 

Maximum  soil  Infiltration 
rate 

Starting  soil  moisture 

Starting  Initial 
abstraction  capacity 

Unit  hydrograph  peak  rate 
factor 

Time  of  concentration 


Technical  literature  supplemented 
by  advice  from  local  SCS 
hydrologists 


Data 


Source 


c.  Snowmelt  (degree-day  method) 

Melt  temperature  threshold 

Starting  snowpack  water 
equivalent 

Melt  rate  coefficient  _ 

Average  dally  or  max  and 
min  dally  temperature 

2.  Land  Use-Related  Data 

Area  of  watershed 

Percent  of  watershed  In  each 
land  use 

♦Percent  Imperviousness  for 
each  land  use 

♦Street  gutter  density 

♦Street  sweeping  Interval 
for  each  land  use 

3.  Mater  Quality-Related  Data 

a.  Surface  Pollutant  Data 

♦Oust  and  dirt  accumulation 
rate  for  each  land  use 

♦Pollutant  fractions  of  the 
dust  and  dirt  (fractions 
of  suspended  solids, 
settleable  solids,  BOD, 
total  nitrogen,  total 
orthophosphate  and  total 
col  1  form 


Hydrologists 

U.S.  National  Weather  Service 

Topographic  maps 
Land  use  maps 

Aerial  photographs  or  field 
surveys 

City  street  maps 

City  Department  of  Public  Works 


Street  sweeping  program 

Lab  analysis  of  the  dust  and  dirt 


Pollutant  accumulation  rates  Field  sampling  program 
b.  Storm  Runoff  Data 


Concentrations  of  above 
pollutants  In  storm  runoff 
for  several  storms 


Field  sampling  program 


Data 


Source 


c.  Dry  Weather  Flow 


Population  is  the  minimum  data 
requirement,  but  several  options 
are  available  for  computing  the 
characteristics  of  dry  weather 

flow.  The  options  require  City  Department  of  Public  Works 

certain  data  on  quantity  and 

quality  of  domestic.  Industrial, 

commercial,  and  pipe 

Infiltration  flows 


4.  Soli  Erosion  (Universal  Soil  Loss 
Equation) 

Soil  classifications  for  all  soils 
in  the  watershed  (Including  slope} 

Soil  erodlbllity  factors  for  each 
soil  type 

Ratio  of  maximum  hourly  rainfall 
intensity  to  the  maximum  30  minute 
intensity 

Overland  flow  erosion  distance 
Ground  cover  factors 
Erosion  control  factors 
Sediment  delivery  ratios 


Soli  maps 

SCS  soil  scientists 

Technical  literature 

SCS  soil  scientists 
SCS  soil  scientists 
SCS  soil  scientists 
SCS  geologists 


* 
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Required  only  for  dust  and  dirt  method  of  pollutant  accumulation 


2.  Discussion  of  Data  Sources  and  Input  Variables. 

a.  Hydrologic  Data.  The  main  block  of  data  required  by  STORM 
Is  the  hourly  precipitation  record.  The  model  can  accomodate  an  unlimited 
amount  of  hourly  precipitation  data.  A  separate  gage  may  be  used  for  each 
subwatershed.  Normally,  ten  to  twenty  years  of  data  will  contain  adequate 
statistical  representation  of  a  given  subwatershed.  Several  gages  may  have 
to  be  used  In  a  study  of  a  large  area  where  rainfall  patterns  are  known  to 
have  spatial  variation.  Since  an  Individual  gage  record  represents  point 
rainfall,  the  rainfall  record  may  have  to  be  corrected  so  as  to  represent 
basin  average  precipitation.  A  single  coefficient  Is  available  In  the 
program  for  correcting  the  rain  record.  If  the  single  correction  coefficient 
Is  not  judged  to  be  adequate,  ihe  user  can  develop  his  own  basin  average 
precipitation  record  to  be  Input  to  the  model.  Precipitation  data  are 
available  on  magnetic  tape,  from  the  U.S.  National  Weather  Service  (Format 
Ho.  488)  for  stations  In  the  U.S.  Otherwise,  the  hourly  data  must  be 
punched  on  cards  in  the  format  for  STORM. 

There  will  often  be  situations  that  existed  which  caused  the  rain 
gage  to  be  Inoperative  for  a  few  hours  to  a  few  days.  This  condition  Is 
"flagged"  by  a  special  character  on  the  tape  In  the  hours  that  It  occurred. 

In  addition,  a  value  will  be  found  at  the  end  of  the  sequence  of  "flagged" 
hours  that  represents  the  total  precipitation  that  occurred  during  the  gage 
failure.  These  situations  must  be  corrected.  Not  only  would  It  be  Incorrect 
to  allow  a  large  amount  of  precipitation  to  remain  in  one  hour,  but  also  the 
special  characters  would  not  be  accepted  by  the  program  format  under  which 
the  rainfall  data  is  read.  The  rain  data  must  also  be  placed  In  the  format 
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described  In  the  users  manual.  The  HEC  has  developed  a  preprocessor 
program  wnich  will  read  the  original  tape,  correct  the  gage  failures  by 
distributing  the  precipitation  evenly  among  the  hours,  reformat  the  data, 
and  place  the  data  In  a  disc-storage  location.  This  allows  It  to  be  called 
up  quickly  without  having  to  go  back  to  the  original  tape. 

uepresslon  storage  (initial  abstraction)  Is  expressed  in  terms  of 
incnes  of  water  depth  over  the  entire  watershed  for  the  coefficient  method 
(each  land  use  for  the  SCS  method).  It  Is  assumed  to  represent  the  sum  of 
the  amount  of  storage  available  In  small  depressions  in  the  land  surface, 
the  amount  Intercepted  by  vegetative  cover,  and  the  amount  of  Infiltration 
that  occurs  during  the  time  when  depression  storage  is  being  filled. 

The  runoff  coefficients  In  the  coefficient  method  are  very  Important 
parameters  since  they  govern  the  quantity  of  runoff  (and  the  quantity 
governs  the  quality  of  runoff).  The  values  should  represent  average  yearly 
fractions  of  rainfall  (in  excess  of  depression  storage)  that  runs  off  the 
watershed.  In  the  coefficient  method,  the  Impervious  area  runoff  rate  Is 
used  to  “drive"  the  pollutant  washoff  equations. 

The  variables  required  for  the  SCS  runoff  method  are  used  for  continuous 
soil  moisture  accounting.  All  are  Important  since  they  govern  the  manner 
In  which  precipitation  losses  are  computed.  Assistance  In  selection  of  the 
magnitude  of  these  variables  for  a  given  watershed  can  be  found  In  References 
3  and  4.  The  most  important  variables  (SMAX,  DEPR,  RATEIH,  and  PERCMX)  can 
normally  be  obtained  from  SCS  soil  survey  reports. 

The  snowmelt  option  requires  four  easily  obtainable  variables:  air 
temperature,  melt  temperature,  melt  rate  coefficient  and  starting  snow  pack 
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water  equivalent.  The  program  compares  each  day's  average  temperature  with 
the  temperature  threshold  In  order  to  decide  whether  precipitation  falling 
on  that  day  is  rain  or  snow.  If  It  Is  snow.  It  goes  Into  the  pack  where  It 
remains  until  a  day  occurs  having  a  mean  temperature  above  the  threshold. 

Melt  is  then  computed  by  applying  the  melt  rate  coefficient  to  the  difference 
between  the  mean  temperature  and  the  threshold  temperature. 

b.  Land  Use-Related  Data.  The  land  use-related  data  make  up  the 
second  main  category  of  Input.  The  land  use  data  are  usually  available  from 
local  planning  agencies.  In  the  San  Francisco  study  (Reference  5),  a 
computerized  system  was  implemented  whereby  the  land  use  for  each  subbasin 
In  the  study  is  produced  by  a  program  which  searches  a  magnetic  tape  of 
digitized  land  use  Information  gathered  by  a  satellite  (LANDSAT).  This 
technique  should  greatly  simplify  the  process  of  obtaining  existing  land 
use  for  a  storm  water  runoff  study. 

The  percentage  Impervious  area  Is  an  Important  variable  In  the 
coefficient  method  of  computing  runoff.  A  distinction  Is  to  be  made  between 
gross  impervious  area  and  connected  Impervious  area.  The  latter  Is  the 
variable  that  should  be  used  In  the  model  since  It  Is  most  highly  correlated 
with  Immediate  runoff  quantity  and  quality.  An  example  of  connected 
Impervious  areas  would  be  direct  connection  of  a  rooftop,  driveway  and 
gutter  to  a  storm  sewer  pipe  system.  Impervious  areas  which  drain  onto 
pervious  areas  should  not  be  included  In  this  variable.  Instead,  the 
runoff  coefficient  for  the  pervious  areas  should  be  Increased  to  account 
for  the  additional  runoff  due  to  non-connected  Impervious  area.  If  gross 
Impervious  area  Is  available  from  techniques  which  use  remotely  sensed 


Information,  It  would  be  desirable  to  develop  a  correlation  between 
connected  and  gross  Impervious  area  by  field  studies. 

c.  Mater  Quality-Related  Data.  The  surface  pollutant  data  required 
are  dependent  on  the  option  chosen  for  pollutant  accumulation.  The  first 
option  (IPACUM*1)  Is  Intended  for  use  In  an  urbanized  watershed  where  the 
majority  of  pollutants  can  be  assumed  to  come  from  street  accumulation  of 
dust  and  dirt.  This  option  requires  a  dust  and  dirt  accumulation  rate  for 
each  land  use,  as  well  as  pollutant  fractions  of  the  dust  and  dirt.  The 
second  option  (IPACUM-2)  requires  only  a  pollutant  accumulation  rate  for 
each  land  use.  It  Is  Intended  to  be  used  In  areas  where  a  significant 
portion  of  the  pollutant  discharge  comes  from  sources  other  than  streets. 

An  example  would  be  a  basin  composed  of  both  urban  and  nonurban  land  uses 
or  a  totally  nonurban  basin.  This  option  could  also  be  used  on  an  entirely 
urban  area  since  the  dust  and  dirt  accumulation  rate  and  pollutant  fractions 
can  be  converted  to  pollutant  accumulation  rates. 

The  second  option  for  surface  pollutant  accumulation  does  not  require 
the  collection  of  street  contaminant  data.  This  option  relies  entirely  on 
using  storm  runoff  data  to  compute  annual  pollutant  discharges  which  are 
used  to  obtain  short  term  accumulation  rates.  The  advantages  of  this  option 
are  that  It  does  not  require  street  contaminant  data  and  It  can  be  used 
where  the  street  surfaces  are  not  believed  to  be  the  major  contributor  of 
pollutants.  Its  disadvantage  Is  that  a  larger  amount  of  storm  runoff  data 
must  be  collected. 

Storm  runoff  data  are  Important  for  calibration  of  the  model.  The 
data  are  required  to  estimate  the  pollutograph  for  Individual  runoff  events. 
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as  well  as  the  annual  discharge  of  a  given  pollutant.  This  topic  will  be 
covered  In  detail  in  Chapter  VI. 

Four  options  are  available  for  computing  the  quantity  and  quality  of 
dry  weather  flow.  The  options  range  from  simply  specifying  the  total 
population  for  the  watershed  to  inputting  the  hourly  variation  of  the  quantity 
(as  well  as  quality)  of  the  dry  weather  flow. 

d.  Soil  Erosion  Data.  Input  data  for  the  land  surface  erosion  comes 
from  several  sources.  The  soil  type  is  normally  available  from  soil  maps. 

The  program  will  allow  use  of  a  representative  sample  of  watershed  area  to 
avoid  repetitious  manual  coding  of  soil  type  Information.  The  ratio  of  30 
minute  to  hourly  rainfall  intensity  is  available  from  technical  literature 
or  the  National  Weather  Service.  The  other  variables  for  this  option  are 
rather  specialized  in  nature.  Therefore,  the  local  Soil  Conservation  Service 
office  should  be  consulted  to  provide  assistance  In  selecting  the  magnitude 
of  these  variables  as  well  as  assistance  in  Interpretation  of  results  from 
the  soil  erosion  analysis. 

The  land  surface  erosion  option  Is  to  be  used  only  for  sediment 
production  studies.  Sediment  loads  calculated  by  this  option  are  not  added 
to  the  suspended  or  settleable  solids  loads  and  therefore  are  not  reflected 
in  the  pollutograph,  event  or  annual  values  of  suspended  and  settleable 
solids.  In  studies  where  soil  erosion  m^y  be  a  contributor,  but  not 
necessarily  the  major  source,  the  loading  coefficients  for  suspended  and 
settleable  solids  must  be  adjusted  In  order  for  the  soil  erosion  to  be 
reflected  in  the  quality  output. 
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VI .  COLLECTION  AND  MANAGEMENT  OF  FIELD  DATA  REQUIRED  FOR  CALIBRATION 

1.  Pollutant  Data.  The  two  types  of  data  required  are  average  annual 
dust  and  dirt  accumulation  rates  and  the  pollutant  composition  of  the  dust 
and  ulrt.  Data  should  be  taken  from  representative  portions  of  each  land 
use.  A  representative  length  (perhaps  100  feet)  of  guttered  street  should 
be  selected  for  study  In  each  land  use.  Criteria  to  judge  Its  representative¬ 
ness  could  Include  average  traffic  volume,  average  population  (or  housing) 
density,  average  age  of  structures,  or  others  considered  to  be  relevant. 

The  measurements  should  be  made  during  a  dry  season.  A  dry  day  accumulation 
Interval  should  be  selected  (perhaps  2  weeks).  The  entire  length  (and 
width)  of  the  100  foot  section  of  street  should  be  "swept".  A  vacuum  device, 
such  as  that  used  by  shopping  centers  to  clean  parking  lots,  may  be  used. 
Brooming  or  washing  may  have  to  be  employed  to  Insure  that  all  fine  material 
is  removed. 

The  first  sweeping  at  each  site  cannot  be  used  to  calculate  accumulation 
rates  since  the  accumulation  time  Is  unknown,  although  It  can  be  analyzed  for 
pollutant  fractions.  Significant  runoff  cannot  be  tolerated  during  the 
accumulation  intervals.  The  field  team  will  have  to  be  alert  to  the 
possibility  that  runoff  could  occur  during  any  accumulation  Interval  and 
thus  invalidate  the  results  for  that  particular  Interval.  The  procedure 
must  then  start  anew,  discarding  the  first  sweeping  and  then  sweeping  again 
at  the  end  of  the  accumulation  period. 

The  material  collected  after  each  accumulation  period  should  be 
screened  over  a  1/4  Inch  mesh  screen.  The  material  passing  the  screen 
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should  be  weighed.  A  sample  of  the  material  should  be  sent  to  a  laboratory 
for  analysis  of  moisture  content,  fractions  of  suspended  solids,  settleable 
solids,  BOD,  total  nitrogen,  total  orthophosphate  or  other  pollutants.  At 
least  6  sweepings  should  be  accomplished  during  the  dry  season.  The  field 
operations  should  be  coordinated  with  city  officials  so  as  to  obtain  approval 
and/or  assistance,  such  as  providing  traffic  barricades  for  use  during 
sweeping.  The  above  data  are  required  for  IPACUM-1. 

2.  Storm  Runoff  Data.  The  purpose  of  this  section  Is  to  outline 
a  procedure  for  collecting  storm  runoff  data  to  calibrate  STORM  on  a  given 
area.  The  basic  premise  Is  to  collect  sufficient  data  to  adequately  define 
the  pollutographs  for  the  required  pollutants  for  enough  storms  to  estimate 
the  annual  pollutant  discharge.  Figure  2  shows  concentrations  of  two 
pollutants  measured  during  a  storm  runoff  event.  Figure  3  shows  a  polluto- 
graph  for  BOD  for  a  single  event. 

In  order  to  obtain  a  better  definition  of  the  quality  of  runoff, 
small  portions  of  the  study  area  should  be  selected  for  detailed  analysis.  In 
this  approach,  a  test  watershed  Is  chosen  to  represent  each  of  the  major 
land  uses  In  the  study  area.  The  model  can  then  be  calibrated  for  each 
land  use,  thereby  evaluating  those  coefficients  which  are  a  function  of 
land  use.  A  major  assumption  of  this  approach  is  that  results  from  a  test 
watershed  will  be  transferrable  to  other  areas  of  the  same  land  use.  The 
use  of  test  watersheds  made  up  of  more  than  one  land  use  Is  not  recommended 
because  one  must  make  adjustments  to  the  coefficients  without  the  benefit 
of  knowing  how  they  vary  with  land  use. 
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Test  watersheds  must  be  chosen  with  good  judgement.  It  Is  Important 
to  choose  basins  that  are  small  enough  to  be  reasonably  homogeneous  In 
terms  of  hydrologic,  pollutant  and  land  use  factors,  and  yet  not  so  small 
that  their  response  will  be  too  short  to  allow  collection  of  adequate  data. 
Depending  on  the  time  of  concentration,  basin  sizes  of  5-10  square  miles 
should  be  adequate. 

Basin  average  rainfall  data  are  required  to  operate  STORM.  A  rainfall 
recorder  Installed  within  the  basin  Is  most  desirable.  More  than  one  gage 
may  be  required  to  measure  the  spatial  variation  of  rainfall  across  the 
basin.  If  more  than  one  gage  Is  used,  the  data  from  each  applicable  gage 
can  be  used  to  compute  basin  average  precipitation  (by  methods  such  as  the 
isohyetal  or  Thiessen) since  STORM  accepts  data  from  only  one  rain  gage  per 
subbasin.  The  recorder  should  be  capable  of  recording  as  low  as  10  minute 
amounts  of  precipitation  to  allow  short  Interval  analysis  of  selected 
events.  It  would  be  desirable  to  have  several  years  of  hourly  rainfall 
data.  The  rain  gage  should  be  installed  according  to  procedures  in 
Reference  6. 

Runoff  quantity  data  are  required  to  calibrate  the  runoff  prediction 
in  the  model.  Mater  quality  predictions  also  depend  on  adequate  repre¬ 
sentation  of  the  runoff  process.  A  streamflow  stage  recorder  must  be 
installed  at  the  outlet  of  each  test  watershed.  Normally,  the  nearest 
U.S.  Geological  Survey  office  can  be  employed  to  Install  and  operate  the 
gaging  stations.  A  reliable  relationship  between  flow  rate  and  stage 
should  be  developed  prior  to  the  study.  Rating  the  gage  may  take  several 
months  depending  on  the  runoff  characteristics  of  the  test  watersheds. 
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If  the  U.S.  Geological  Survey  cannot  Install  and  operate  the  gaging  stations. 
It  should  be  accomplished  according  to  procedures  outlined  In  Reference  7. 

Adequate  storm  water  quality  data  are  the  most  difficult  and  costly 
to  obtain.  Ideally,  all  significant  runoff  events  during  (me  or  more  years 
should  be  measured  so  as  to  obtain  estimates  of  the  average  annual  discharge 
of  selected  constituents.  Sampling  costs  usually  govern  the  number  of 
storms  that  can  be  measured,  therefore  the  researcher  must  be  satisfied 
with  measurements  from  a  small  number  of  storms.  In  this  case.  It  may  be 
possible  to  develop  a  relationship  between  pollutant  discharge  and  runoff 
Intensity,  number  of  dry  days  prior  to  the  event,  or  other  relevant  variables. 
This  procedure  can  then  be  used  to  estimate  the  average  annual  loading  of  a 
given  pollutant.  In  the  San  Francisco  study  (Reference  5),  It  was  not 
possible  to  choose  which  storms  were  to  be  monitored,  a  priori.  Due  to 
the  extremely  erratic  nature  of  runoff  occurrences  In  that  area,  the  re¬ 
searchers  had  to  attempt  to  measure  each  runoff  event  In  hopes  that  It 
would  be  significant. 

Sampling  techniques  range  from  completely  manual  to  completely 
automatic,  bach  technique  has  met  with  limited  success  depending  on  the 
procedure,  equipment  and  rainfall -runoff  characteristics  of  the  test  water¬ 
sheds.  In  the  manual  approach,  the  field  team  must  rely  upon  weather 
forecasts  to  prepare  for  sampling  an  event.  At  the  onset  of  precipitation, 
tne  team  must  be  imnediately  dispatched  to  the  sampling  site.  Sampling  must 
begin  imnediately  so  as  to  obtain  measurements  that  represent  the  'first 
flush"  or  Initial  washoff  phenomenon.  Sampling  must  continue  at  frequent 
Intervals  (say  15  minutes)  until  the  peak  Is  reached  since  the  majority 
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of  the  mass  of  pollutants  are  assumed  to  be  washed  off  during  this  period. 
When  the  hydrograph  peak  Is  reached  (this  may  be  difficult  to  recognize 
since  many  storms  produce  multi-peaked  hydrographs)  the  sampling  frequency 
can  be  reduced  (to  say,  1  hour). 

The  opposite  extreme  Is  using  a  fully  automated  approach.  Equipment 
is  available  to  pump  many  samples  from  the  stream  Into  separate  sample 
bottles  at  specified  Intervals,  usually  5-60  minutes.  Some  devices,  such 
as  that  used  by  Florida  Institute  of  Technology,  (Reference  8)  can 
refrigerate  the  samples  for  up  to  48  hours  to  preserve  the  samples  before 
they  are  retrieved  by  the  field  crew.  These  devices  can  be  programmed  to 
begin  sampling  at  a  significant  rise  in  stage,  and  some  will  sample  at 
preset  increments  of  either  time  or  stage.  Most  automatic  sampling 
equipment  that  has  been  used  Is  of  the  type  that  conmences  and  samples 
automatically,  but  relies  upon  a  field  crew  to  retrieve  the  samples  within 
l-<!  hours. 

The  automatic  sampling  technique  has  had  varied  success  In  storm 
water  applications.  The  equipment  seems  to  be  prone  to  mechanical  failure. 
Runoff  is  Intermittent  requiring  that  the  equipment  lie  Idle  for  some  time 
and  yet  function  properly  at  the  beginning  of  a  runoff  event.  The  situation 
is  compounded  by  the  nature  of  storm  water.  The  flow  rate  Is  extremely 
variable,  and  the  water  is  often  laden  with  debris  which  causes  clogging 
of  suction  lines  or  pumps.  In  the  San  Francisco  Study  (Reference  5), 
automatic  sampling  at  the  Castro  Valley  watershed  was  discontinued  In  1973 
because  of  mechanical  failure  and  clogging.  Since  then,  the  field  team 
has  reliea  upon  reaching  the  sampling  site  within  30  minutes  after  runoff 


begins.  Part  of  the  Initial  runoff  was  sacrificed  In  order  to  Insure  that 
samples  were  collected  throughout  the  runoff  event.  However,  In  the  Strong 
Ranch  watershed,  the  automatic  equipment  was  retained  since  the  site  Is 
within  walking  distance  of  the  field  office,  and  a  technician  was  able  to 
Insure  that  the  equipment  would  begin  sampling  at  the  correct  time  and 
continue  sampling  throughout  the  event. 

Sampling  must  be  accomplished  so  as  to  provide  enough  data  points 
to  adequately  define  the  pollutographs  for  each  selected  constituent. 

The  pollutographs  are  Integrated  to  determine  the  total  discharge  of 
pollutants  by  each  measured  event.  Sampling  frequency  and  total  duration 
of  sampling  are  Important  factors  In  defining  the  pollutographs.  A  rule 
of  thumb  Is  that  4-6  points  should  be  defined  on  the  rising  limb  of  the 
hydrograph  and  2-4  on  the  descending  limb.  If  6  points  on  the  rising  limb 
are  chosen,  the  sampling  frequency  should  be  approximately  1/6  the  average 
time  to  peak  for  the  basin.  When  using  automatic  equipment,  more  samples 
than  this  can  easily  be  collected.  Figure  3  shows  a  good  distribution  of 
sampling  points  throughout  a  pollutograph. 

Regardless  of  the  sampling  method,  the  samples  must  be  handled 
properly.  The  sample  bottles  should  be  placed  In  a  mixture  of  Ice  and  water 
as  soon  as  possible  after  collection.  This  Is  particularly  Important  for 
samples  that  will  be  analyzed  for  biological  parameters  such  as  BOD.  These 
samples  should  be  analyzed  within  24  hours.  Most  physical  and  chemical 
constituents  can  be  analyzed  later.  Laboratory  analysis  should  be  performed 
according  to  Standard  Methods  for  Examination  of  Wastewater  (Reference  9). 
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Since  the  present  version  of  STORM  predicts  6  parameters  (suspended 
solids,  settleable  solids,  5-day  BOO,  total  nitrogen,  total  orthophosphate, 
and  total  conform  bacteria)  It  Is  Imperative  that  tests  be  made  for  these. 
Certain  other  parameters  such  as  pH,  alkalinity,  dissolved  oxygen  and  total 
dissolved  solids  should  also  be  measured  regularly.  Other  parameters  such 
as  heavy  metals,  COO,  oil  and  grease,  fecal  conforms,  herbicides  and 
pesticides  can  also  be  measured  If  time  and  funds  permit.  It  may  be 
possible  to  predict  these  additional  parameters  by  adjusting  the  coefficients 
and  equations  In  STORM  used  for  the  other  constituents.  If  the  parameters 
cannot  be  related  to  the  dust  and  dirt,  sufficient  data  must  be  collected 
to  establish  a  functional  relationship  between  the  washoff  of  that  pollutant 
and  measurable  parameters  such  as  runoff  Intensity,  number  of  dry  days 
between  events,  traffic  volume  or  other  Independent  variables.  The  maximum 
number  of  pollutants  predicted  during  any  run  cannot  exceed  6. 

Efficient  data  handling  can  greatly  simplify  a  storm  water  study. 

Since  data  are  to  be  taken  from  3  main  sources,  (rainfall,  runoff  quantity 
ana  runoff  quality)  It  Is  Imperative  that  data  coordination  be  handled  by 
a  single  Informed  entity.  For  example,  pollutant  concentrations  are  not 
very  meaningful  In  themselves,  but  must  be  viewed  together  with  flow  rate 
and  time.  Data  can  most  efficiently  be  analyzed  by  placing  them  on  cards 
so  that  analysis  and  plotting  of  data  and  results  can  be  handled  by 
computer.  The  HEC  has  developed  a  format  for  these  cards  which  Include 
precipitation,  flow,  and  concentration  data.  The  format  Is  shown  In 
Figure  4.  An  example  of  output  from  a  pollutograph  computation  and  plotting 
routine  developed  by  HEC  Is  shown  In  Figure  5. 
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The  data  can  also  be  submitted  to  the  Environmental  Protection  Agency 
for  placing  In  their  STORET  system  for  storage  and  retrieval  of  water  quality 
aata.  This  makes  the  data  available  to  anyone  In  the  U.S.  who  has  access  to 
the  system.  Certain  plotting  and  statistical  analysis  capabilities  are 
available  during  an  Interrogation  of  the  system.  The  procedures  for  accessing 
the  STORET  system  are  outlined  In  Reference  10.  An  example  of  data  from  a 
STORET  retrieval  Is  shown  In  Table  1. 

VII.  MODEL  CALIBRATION 

1.  Introducti on.  Calibration  can  be  characterized  as  a  "mini 
application".  Since  the  model  must  be  applied  to  the  calibration  watersheds, 
all  those  variables  required  for  a  normal  application  must  be  assembled  for 
use  In  calibration.  The  procedure  entails  assemblage  of  observed  data 
(rainfall,  runoff  and  quality)  and  land  use  data  followed  by  successive 
applications  of  the  model  to  the  test  watersheds  until  the  model  Is 
satisfactorily  calibrated  for  the  study  area. 

2.  Uuanti ty .  Model  calibration  cannot  commence  until  certain  data 
are  prepared  for  input  to  the  model  and  observed  data  organized  so  as  to 
compare  with  computed  results. 

Precipitation  data  must  be  extracted  from  recorder  charts.  These 
charts  usually  contain  traces  which  represent  accumulated  precipitation. 

The  value  at  the  end  of  each  hour  can  be  tabulated.  The  Incremental  amounts 
are  found  by  subtracting  the  preceding  value  from  each  hour's  accumulated 
total.  The  hourly  precipitation  must  then  be  coded  In  the  format  described 
In  the  users  manual  (Reference  2). 
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The  next  step  Is  assembling  the  runoff  quantity  data.  If  the  U.S. 
Geological  Survey  handles  the  streamflow  gaging,  the  data  will  normally  be 
furnished  In  the  form  of  mean  dally  discharges  plus  selected  Individual 
storm  hydrographs.  Monthly  and  annual  runoff  amounts  are  also  furnished. 
Because  processing  and  publishing  the  flow  of  data  often  takes  several 
months,  sufficient  time  must  be  allowed  for  this  Item.  If  the  stream  gaging 
is  handled  "In-house",  the  stage  recorder  charts  must  be  converted  to  flow 
rates  by  use  of  the  rating  curve  (previously  established)  for  each  Individual 
gaging  station.  The  hydrographs  are  then  Integrated  so  as  to  produce 
statistics  such  as  individual  hydrograph  runoff  volumes,  monthly  runoff 
volumes,  and  annual  volumes. 

The  quantity  calibration  Involves  adjustment  of  certain  coefficients 
which  regulate  the  volume  and  timing  of  runoff.  In  the  coefficient  method, 
the  first  factor  to  adjust  Is  the  runoff  coefficient  for  the  pervious  areas 
of  the  watershed  (CPERV)  so  that  the  observed  and  computed  magnitudes  of 
total  period  runoff  volume,  annual  volume  and  monthly  volume  show  fair 
agreement.  Since  the  coefficient  method  Is  a  very  crude  model  of  the  urban 
runoff  process,  it  may  not  be  meaningful  to  compare  individual  event  volumes. 
Other  variables  that  can  be  adjusted,  although  to  a  lesser  extent,  are  the 
runoff  coefficient  for  the  Impervious  areas  (CIMP)  and  the  initial  abstraction 
(OEPRS). 

The  SCS  runoff  method  should  produce  more  accurate  magnitudes  of 
runoff  excesses.  Therefore,  when  using  this  method,  one  should  not  only 
produce  more  accurate  total  period,  annual,  and  monthly  volumes  of  runoff, 
but  also  more  accurate  Individual  hydrograph  volumes. 
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The  unit  hydrograph  procedure  should  Increase  the  accuracy  of  the 
timing  and  shape  of  the  Individual  hydrographs.  Table  2  shows  the  result 
of  preliminary  sensitivity  tests  using  the  SCS  methods. 

3.  quality.  Runoff  quality  data  must  also  be  assembled.  A  tabular 
listing  should  be  made  of  the  date,  time  and  concentrations  In  mg/1  for  the 
6  parameters  predicted  by  STORM  (conforms  in  1000  MPN/1).  Table  1  shows 
sucn  a  listing.  Pollutographs  for  each  parameter  should  be  computed  by 
multiplying  the  concentrations  by  the  water  discharge  rate  at  the  time  of 
measurement.  The  values  should  be  converted  to  Ibs/hr.  The  pollutographs 
should  be  integrated  so  as  to  obtain  the  total  pollutant  discharge  for  each 
runoff  event. 

Results  from  the  street  sweeping  program  must  be  expressed  In 
appropriate  units  for  input  to  the  model.  The  dust  and  dirt  collected 
at  the  end  of  each  accumulation  Interval  should  be  divided  by  the  Interval, 
averaged,  and  expressed  on  a  dry  weight  basis  for  each  land  use  (lbs/100  ft. 
gutter/day).  The  constituents  of  the  dust  and  dirt  should  be  expressed  In 
terms  of  average  pounds  of  pollutant  per  100  pounds  of  dust  and  dirt  for 
each  land  use.  This  will  provide  the  data  for  use  of  pollutant  accumulation 
method  number  1  (IPACUM=1). 

If  street  sweeping  data  are  not  available  or  If  the  streets  are  not 
believed  to  be  the  major  source  of  pollutants,  the  second  method  for 
pollutant  accumulation  must  be  used  (IPACUM»2).  This  method  relies  upon 
the  stormwater  runoff  data  alone  to  establish  the  accumulation  rate  for 
each  pollutant.  Several  assumptions  must  be  made  In  order  to  express  the 
pollutant  accumulation  rates  in  terms  of  pounds/acre/day.  The  first,  and 


31 


observed 


computed 


L  ■  little  effect,  M  ■  moderate  effect, 

H  ■  high  effect 

I  ■  Increases  runoff,  D  *  decreases  runoff 
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probably  the  most  significant.  Is  that  the  annual  mass  discharge  of  a  given 
pollutant  can  be  estimated  by  measurement  of  the  mass  discharges  for  several 
selected  events  and  extrapolation  to  an  entire  year.  The  extrapolated  can 
be  done  by  correlation  with  relevant  washoff  variables  such  as  water  runoff 
volume  and  rates.  The  second  assumption  Is  that  the  dally  accumulation  rate 
can  be  estimated  by  dividing  the  total  annual  mass  discharge  by  the  average 
number  of  dry  days  per  year. 

It  Is  recognized  that  the  accumulation  rates  must  vary  seasonably. 
However,  no  data  are  available  to  quantify  this  variability.  Therefore, 
one  must  proceed  with  a  single  linear  dally  accumulation  rate.  A  third 
assumption  Is  that  the  net  change  In  annual  mass  discharge  of  pollutants 
Is  negligible. 

The  state-of-the-art  In  urban  storm  water  quality  modeling  precludes 
highly  accurate  simulation  of  pollutographs.  However,  the  quality  aspects 
of  STORM  must  be  calibrated  so  as  to  produce  results  as  accurate  as  possible. 
It  may  not  be  possible  to  reproduce  the  time  value  of  concentrations  In 
which  case  the  emphasis  should  be  placed  on  reproducing  the  mean  concentration 
and  the  total  mass  of  pollutants  washed  off  by  the  observed  events. 

Quality  calibration  can  commence  only  after  an  adequate  quantity 
calibration  has  been  completed.  An  initial  trial  can  be  made  using  the 
suggested  values  for  pollutant  accumulation  (lbs/100  ft  gutter/day  or 
Ibs/acre/day).  The  event  numbers  to  be  studied  In  detail  can  be  taken 
from  this  Initial  run  and  a  preliminary  Inspection  of  the  extent  of 
agreement  between  observed  and  computed  data  can  be  made.  Comparisons 
should  be  made  first  with  observed  and  computed  pollutant  concentrations 


(displayed  In  the  pollutograph  analysis  for  the  selected  observed  events). 

The  primary  coefficients  to  be  adjusted  are  the  pollutant  accunulatlon 
rates  and  the  washoff  exponent,  the  latter  being  a  coefficient  which  governs 
the  rate  of  washoff  of  each  pollutant.  In  the  dust  and  dirt  method,  the 
dust  and  dirt  accumulation  rate  for  each  land  use  can  be  adjusted  first 
depending  on  whether  all  pollutant  concentrations  are  too  low  or  too  high. 
Once  the  concentrations  are  nearly  correct,  the  pollutant  factors  of  the 
dust  and  dirt  can  be  adjusted.  (Ideally,  the  model  should  be  calibrated 
on  single  land  use  basins,  as  recommended  earlier.)  The  washoff  exponent 
can  also  be  adjusted  If  the  data  show  a  different  change  In  concentration 
with  time  during  each  event.  (A  higher  exponent  will  produce  a  higher  rate 
of  washoff.)  In  the  dally  pollutant  accumulation  method,  the  accumulation 
rates  for  each  pollutant  and  the  washoff  exponent  are  adjusted  until  fair 
agreement  exists  between  the  measured  and  computed  data.  After  agreement 
in  concentrations  is  accomplished,  one  should  compare  the  total  mass  of 
each  pollutant  washed  off  with  the  computed  values  for  the  selected  events. 
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IX.  MODEL  APPLICATION 


1.  General .  This  section  presumes  that  STORM  has  been  properly 
calibrated  on  the  test  watersheds  (both  quantity  and  quality).  The 
remaining  tasks  Involve  application  to  the  entire  study  area. 

The  first  step  In  an  application  Is  to  delineate  the  entire  study 
area  on  a  map.  A  good  scale  to  use  Is  the  1:24,000  size  of  the  U.S. 
Geological  Survey.  The  study  should  be  outlined  by  drainage  areas.  The 
study  area  must  be  subdivided  (by  drainage,  again)  for  application  of  the 
model.  The  subdivision  is  required  to  generate  loadings  at  correct  locations 
along  the  receiving  stream,  and  for  routing  flows.  A  usable  number  of 
subbasins  is  25-30  since  a  data  deck  must  be  prepared  to  each  subbasin. 

Since  the  program  has  the  capability  to  do  several  subbasins  In  one  run, 
it  is  recommended  that  the  receiving  stream  be  subdivided  into  a  number  of 
subreaches,  say,  4-5.  The  loading  points  for  each  reach  can  be  run  at  the 
same  time.  If  computational  or  storage  capability  is  not  available  to  run 
several  subbasins  at  a  time,  they  may  be  Individually  run.  In  either  case 
the  STORM  output  is  stored  on  tape/disk  for  later  Input  to  a  receiving 
water  assessment  model. 

If  rainfall  gage  locations  other  than  those  used  in  the  calibration 
phase  are  to  be  used,  that  data  must  be  assembled  for  use  in  the  application. 
The  rain  data  can  usually  be  obtained  on  magnetic  tape  from  the  National 
Weather  Service  in  Asheville,  North  Carolina.  Some  editing  of  the  data 
supplied  on  tape  may  be  required,  as  discussed  In  Chapter  V.  Once  the 
point  rainfall  data  have  been  satisfactorily  edited,  reformatted,  and 
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placed  In  a  usable  storage  location,  a  decision  must  be  made  In  regard  to 
which  gages  to  associate  with  the  various  subbasins.  The  proximity  of  the 
centroid  of  the  basin  to  the  rain  station  Is  probably  the  major  consideration. 
A  decision  must  also  be  made  as  to  whether  to  use  the  point  rainfall  or  to 
estimate  subbasin  average  precipitation.  STORM  has  the  capability  of 
multiplying  each  hour's  precipitation  amount  by  a  constant  (specified  by 
the  user)  for  use  In  situations  where  this  Is  deemed  adequate. 

The  next  major  block  of  required  Input  data  Is  the  land  use  Information. 
Land  uses  are  specified  In  terms  of  the  percentage  of  each  In  the  subbasin. 

In  the  Oconee  River  Study  (Reference  11),  the  land  use  percentages  were 
calculated  from  a  computer-stored  data  file  of  basin  land  use  delineated  on 
a  1.5  acre  grid  cell  basis.  STORM  was  modified  to  retrieve  this  Information 
directly  from  tape/disk.  This  avoided  the  expense  of  time  and  funds  for 
placing  the  land  use  Information  on  cards  for  use  by  STORM. 

2.  Pollutographs.  Assuming  that  the  model  has  been  calibrated  on 
single  land  use  basins,  the  application  becomes  a  matter  of  Inputting 
those  factors  which  are  a  function  of  land  use  for  each  land  use  represented 
In  each  subbasin.  A  separate  set  of  E-T  cards  are  to  be  punched  for  each 
subbasin.  The  subbasins  can  be  Identified  by  number  and  river  mile  along 
the  receiving  stream  where  the  subbasin  Is  assumed  to  empty.  The 
subbasin  cards  can  be  grouped  In  reaches  and  those  reaches  run  In  a  single 
pass. 

In  studies  where  the  effects  of  storm  water  runoff  on  the  receiving 
waters  are  assessed,  It  has  been  conventional  to  study  a  "critical"  period 
of  In-stream  water  quality.  The  criteria  for  selection  of  this  critical 
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period  are  somewhat  Indefinite  but  Include  warm  weather  and  low  flow 
conditions.  The  effects  of  several  historical  storm  water  runoff  events 
on  the  receiving  water  are  usually  studied.  This  procedure  was  followed 
In  the  St.  Louis  (Reference  12),  Atlanta  (Reference  13)  and  Oconee  studies 
(Reference  13). 

The  STORM  model  Is  then  used  to  predict  the  wet  weather  storm  water 
runoff  pollutographs  for  the  events  during  the  critical  period.  (It  will 
normally  be  required  to  use  about  S  years  of  precipitation  data  In  STORM 
prior  to  the  selected  events  so  that  the  effects  of  Initial  conditions  can 
be  neglected.)  The  current  version  of  STORM  will  print  the  hourly  polluto- 
graph  for  each  selected  event  (maximum  s  20)  although  it  gives  a  one  line 
summary  of  each  event  that  occurred  during  the  period  of  simulation.  It 
will  normally  be  convenient  to  modify  the  program  slightly  so  as  to  punch 
the  pollutographs  directly  on  cards  or  place  them  on  tape/disk.  This  allows 
the  pollutographs  to  be  used  directly  by  the  receiving  water  assessment 
model  without  any  Intermediate  data  manipulations  by  the  user.  A  future 
version  of  STURM  will  contain  a  simple  In-stream  water  quality  assessment 
module,  as  well  as  a  module  to  route  and  combine  flows  and  quality  from 
subbasins. 

3.  Storage -Treatment  Rate  Analysis.  STORM  has  a  great  deal  of 
potential  for  providing  Information  which  could  assist  In  the  preliminary 
sizing  of  storage  volumes  and  treatment  rates  to  achieve  desired  control 
of  storm  water  runoff.  If  treatment  of  the  stored  water  is  not  envisioned, 
the  treatment  rate  could  be  equated  to  a  pumping  rate  so  that  the  model 
could  be  used  to  study  a  flood  control  reservoir  with  a  constant  release 
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rate.  Minor  modifications  to  the  model  would  allow  the  release  rate  to  be 
made  a  function  of  the  head  on  the  conduit  since  a  variable  release  rate 
Is  more  advantageous  for  optimum  flood  control  operation.  These  options 
will  be  available  In  a  future  version. 

If  values  for  both  storage  and  treatment  are  specified,  the  model 
will  analyze  the  system's  response  to  the  continuous  record  of  runoff 
hydrographs.  The  runoff  amounts  are  computed  by  simulation  of  the  rainfall- 
runoff  process  as  discussed  earlier.  Part  of  the  purpose  for  a  continuous 
analysis  Involves  the  Interaction  of  successive  runoff  events.  From  a  flood 
control  standpoint,  it  would  be  most  advantageous  to  have  the  reservoir  In 
an  empty  state  at  the  beginning  of  each  runoff  event.  This  is  often  not 
possible  from  a  practical  standpoint  and  because  of  the  random  occurrences 
of  storms.  Obviously,  If  the  reservoir  Is  partly  full  from  a  prior  storm. 

It  will  be  less  effective  In  capturing  all  of  a  subsequent  runoff  event. 

STORM  provides  Inormatlon  on  the  characteristics  of  the  stored 
runoff  such  as  the  time  storage  began,  time  of  spill  (or  maximum  storage), 
hours  to  empty,  hours  between  use  of  storage,  to  name  a  few.  Statistics 
are  also  presented  on  the  utilization  of  storage.  These  Include  the  average 
storage  at  each  hour  of  each  event,  the  average  annual  number  of  hours  each 
hundredth  of  an  inch  of  storage  was  utilized.  These  two  statistics  are 
used  to  compute  a  normalized  storage  utilization  curve  which  depicts  the 
percent  of  time  that  the  water  In  storage  was  less  than  or  equal  to  a 
given  percentage  of  the  maximum  storage. 

Figure  0  shows  an  example  of  the  normalized  storage  utilization 
curve.  It  gives  a  quick  visual  measure  of  the  relative  efficiency  of  the 
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selected  storage  volume  and  treatment  rate.  It  gives  an  estimate  of  the 
response  of  the  storage/treatment  combination  to  the  continuous  sequence 
of  runoff  hydrographs.  Referring  to  Figure  6,  one  can  see  that  the 
combination  represented  by  the  dashed  line  Is  not  an  efficient  combination 
since*  for  example,  80%  of  the  time  only  about  20%  of  the  storage  was 
utilized.  Since  the  reservoir  level  would  be  in  the  lower  1/5  of  the 
storage  80%  of  the  time,  we  would  conclude  that  either  the  treatment  rate 
or  the  maximum  storage  capacity  Is  too  high.  The  decision  to  lower  either 
the  storage  or  the  treatment  rate  must  be  made  In  terms  of  the  economics 
of  both  and  the  criteria  for  overflows  from  storage.  The  solid  line  on 
Figure  6  shows  a  more  efficient  combination  of  storage  and  treatment 
rate. 

The  model  Is  designed  to  compute  statistics  for  a  matrix  of  storages 
and  treatment  rates.  Computational  time  Is  the  only  limiting  factor  on  the 
number  of  combinations  that  should  be  analyzed  In  a  given  run.  The 
computer  run  times  are  also  proportional  to  the  length  of  continuous 
record  to  be  studied,  experience  In  HEC  has  shown  that  computational 
time  Is  extremely  sensitive  to  the  combination  of  storage  and  treatment 
rate.  In  fact,  the  computer  time  Is  directly  proportional  to  the  storage 
volume  and  Inversely  proportional  to  the  treatment  rate.  The  reason  for 
the  sensitivity  Is  that  for  the  larger  storages  and  smaller  treatment  rates, 
there  Is  less  spill  and  the  program  must  keep  account  of  more  water  for  a 
longer  period  of  time.  The  HEC  has  added  an  option  In  the  model  so  that 
the  user  may  select  whether  or  not  he  will  compute  the  ages  of  storage 
since  these  Involve  significant  computer  time  for  the  high  storage-low 
treatment  rate  combinations. 
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Figure  7  shows  a  set  of  quantity  performance  curves  which  were 
plotteu  from  an  analysis  of  a  matrix  of  storage  and  treatment  rate  combinations. 
One  can  see  that  If  a  standard  of  not  more  than  2  Inches  overflow  per  year 
was  adopted,  the  shaded  region  on  Figure  7  represents  the  variety  of  storage 
and  treatment  rates  that  would  satisfy  the  standard.  Figure  8  shows  an 
analagous  set  of  curves  for  quality  (BOD,  In  this  case).  The  shaded  region 
represents  a  variety  of  combinations  of  storages  and  treatment  rates  that 
would  release  not  more  than  5000  pounds  of  BOD  per  year  In  overflows  from 
storage. 

Another  standard  could  be  adopted  where  two  or  more  criteria  could 
be  combined  such  that  both  must  be  met  simultaneously.  For  example,  the 
decision  criteria  could  be: 

(1)  not  more  than  5  overflows  per  year,  and 

(2)  not  more  than  2500  pounds  of  BOD  overflow  per  year. 

A  family  of  curves  relating  storage,  treatment  rate  and  average  annual 
number  of  overflows  can  be  prepared  from  the  output  from  STORM.  These 
curves  are  displayed  in  Figure  9  where  a  region  of  less  than  or  equal  to 
2500  pounds  of  BOO  overflow  per  year  has  been  shaded.  It  can  be  seen  that 
point  A  represents  the  amount  of  storage  and  treatment  rate  which  will 
exactly  satisfy  both  criteria.  Point  A  Is  located  at  a  storage  of  0.39 
Inches  and  a  treatment  rate  of  0.033  Inches  per  hour  and  represents  minimum 
requirements  for  both  storage  volume  and  treatment  rate.  The  entire  region 
above  the  5  overflow  line  and  the  2500  pounds  of  BOD  overflow  line  would 
more  than  satisfy  both  criteria.  A  solution  Is  then  pursued  which  would 
satisfy  the  criteria  at  minimum  cost. 
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Quantity  performance  curves 


Figure  8 


Quality  performance  curves 


SHADED  REGION  INDICATES 
BOD  OVERFLOW  $  2500  LBS/YR 


System  performance  curves 


Another  use  of  the  output  from  STORM  Is  to  analyze  the  percent 
removal  of  pollutants  by  virtue  of  detainment  In  the  storage  reservoir. 

The  current  version  does  not  model  the  treatment  process  Itself  nor  the 
settlement  and  decay  of  pollutants  while  In  storage.  Therefore,  the  analysis 
Is  based  entirely  on  detainment  of  the  runoff  until  the  treatment  plant  can 
accommodate  It.  Figure  10  shows  an  example  of  the  percent  of  BOO  removed 
from  runoff  versus  storage  capacity.  In  this  analysis  It  was  assumed  that 
the  treatment  plant  removed  90  percent  of  all  BOD  In  the  Influent  and  the 
BOO  concentration  In  the  overflow  was  the  same  as  the  Influent.  Therefore, 
the  curve  becomes  assymtotlc  to  90  percent  removal  as  storage  Is  Increased 
and  all  runoff  can  be  treated.  A  future  version  of  STORM  will  contain 
algorithms  for  treatment  and  settlement  of  pollutants  In  storage.  This  will 
allow  a  more  realistic  appraisal  of  the  effect  of  storage  and  treatment  on 
the  quality  and  timing  of  runoff  released  to  the  receiving  waters. 

The  preceding  examples  were  only  a  few  Illustrations  of  how  STORM 
can  be  used  to  study  the  quantity  and  quality  of  storm  water  runoff.  Many 
other  specialized  criteria  and  methods  can  be  utilized  as  required  by  the 
application. 
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